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ABS TRACT 


Goda's (1975) model, describing wave transformation from 
deep water to across the surf zone, 1S compared with a large 
amount of wave data obtained from experiments conducted at 
Torrey Pines Beach, San Diego, California. Goda's model 
Simulates wave breaking by truncating the Rayleigh distribu- 
tion in order to estimate the wave height distributions across 
the surf zone; wave heights are shoaled by applying nonlinear 
theory. Comparisons between the empirical distributions and 
theoretical distributions, and between measured and theoreti- 
cal rms wave heights, are made. It is found that Goda's 
model over-predicts the tails and under-predicts the peaks 
of the empirical distributions, and that the calculated rms 
wave heights are too large compared with measured values. 

The range of breaking, and the coefficients used in the 
breaking criteria by Goda, are modified in order to obtain 
a model which better fits the distribution of observed 
heights, and which matches the model and observed rms wave 
heights. The results are quite good, with error envelope 
for predicted rms wave heights less than 20%. Linear shoal- 
ing theory is applied to the model and found to be as good 


as applying nonlinear theory. 
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fee kODUCTION 


The evaluation of an irregular group of shoaling waves 
as they approach and pass through the breaker zone is a com- 
plex process which requires special measurements and analy- 
Sis considerations. The usual approach to shallow water 
wave transformations is to predict, from a single "represen- 
tative" set of deep water parameters, the wave height, the 
wavelength and the frequency at specific shallow water depths, 
using linear shoaling theory. The primary objection to this 
approach is that a single set of deep water wave parameters 
does not realistically represent the distributional charac- 
teristics of naturally occurring sea surface waves. A 
secondary objection arises from the use of linear transfor- 
Mations which become inadequate when applied through the 
surf zone (Wood, 1974). 

Wave heights in deep water, having Gaussian surface ele- 
Vetrems, are described by the Rayleigh distribution (Longuet- 
Higgins, 1952). Waves propagating toward shore can increase 
in height due to shoaling effects, refraction and wave inter- 
actions, and eventually reach a depth where they start break- 
ing. The energy dissipation due to breaking has been simu- 
lated (Goda, 1975) by truncating the tail of the Rayleigh 
Gusteribution. 

Experiments were conducted at Torrey Pines Beach, San 
Diego, California, during November 1978. Sea surface eleva- 


tions, pressures and velocities were measured at closely 
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spaced locations in a line extending from 10m depth to 

inside the surf zone. This thesis applies Goda's model to 
the measurements in order to examine the shoaling and trans- 
formation of wave heights and their probability density func- 
tions (pdf's) from deep water to breaking and across the 


surf zone to the shoreline. 
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II. THEORETICAL BACKGROUND 


A. RAYLEIGH DISTRIBUTION 

The Rayleigh distribution was shown theoretically by 
Longuet-Higgins (1952) to apply to ceep water wave heights 
on the assumption that the sea waves are a narrow-banded 
Gaussian process. Barber (1950) had earlier presented empiri- 
cal evidence that the Rayleigh distribution agreed with the 
measured distribution of waves. On the assumption that the 
wave height is twice the wave amplitude, the wave height 


probability density is then represented by 


atm) 2 2H/He exp (-H*/H~_ (1) 


Ss 5) 


where Hem is the rms wave height. 


Ss 
Using pressure records in the Gulf of Mexico, Longuet- 
Higgins (1975) observed that the Rayleigh distribution fits 
the observed distribution reasonably well in "fairly deep 
water". He found that there is a slight excess of waves 
with heights near the middle of the range and a deficit at 
the two extremes. Since much of the high-frequency portions 
of the wave records were filtered out by the pressure trans- 
ducer, Longuet-Higgins (1975) suggested that the narrow band 
approximation may not be as applicable for the unfiltered 
records. In shallow water with much steeper waves, the 


Rayleigh distribution can again be expected to be less 


applicable due to the non-linearities. 
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Since the Rayleigh distribution theoretically did not 
apply to broadband wave spectra, Goda (1970) numerically 
Simulated wave profiles, where the amplitudes were specified 
by various theoretical spectra of varying bandwidth and the 
phase was random. He then examined the simulated records for 
surface elevations, crest-to-trough wave heights and zero- 
up-crossing wave heights. He found that, using the zero-up- 
crossing determination of wave heights, the Rayleigh distri- 
bution is a good approximation irrespective of the spectral 
bandwidth. Tayfun (1977), in studying the transformation of 
deep water waves to shallow water waves, showed that the 
Rayleigh distribution for wave amplitude was generally 
applicable to all bandwidths. 

The Rayleigh distribution is applied correctly only to 
low waves in deep water (Longuet-Higgins, 1975), since it is 
assumed that the contributions from different parts of the 
generating area are linearly superposable. Under this assumpt- 
tion, the distribution clearly should not hold for waves 
approaching maximum height, i.e., close to breaking, as in 
the surf region or even in the open sea with whitecaps. It 
has been found by several authors (Chakrabarti and Cooley, 
1977; Forristal, 1978) that the theoretical Rayleigh dis- 
tribution over-predicts the maximum wave in the tail compared 
with large wave observations. Forristal (1978) attributed 
the differences to the non-linear, non-Gaussian and skewed 


nature of the free surface. 
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Tayfun (1980) examined non-linear effects by consider- 
ing an amplitude-modulated Stokian wave process with the 
restriction that the underlaying first order spectrum is 
narrow band. The surface displacements were found to be 
non-Gaussian and skewed, and wave heights distributed ac- 
cording to the Rayleigh probability law, particularly for low 
and medium wave height ranges. On the basis of the results 
obtained, Tayfun (1980) concludes that the non-Gaussian 
characteristics of the free surface do not directly result 
in reducing maximum wave heights in a manner consistent with 
field observations and that a more plausible mechanism is 
wave breaking, which is a non-linear effect not directly 
accounted for in the analytical wave models currently available. 

Longuet-Higgins (1980) analyzed the effects of non- 
linearity and finite bandwidth on the distribution of wave 
heights to explain the differences with observations found 
by Forristal. He found that the reason for the discrepancy 
could be accounted for by the presence of free background 
"noise" in the spectrum, outside the dominant peak, and that 
it was not due to a finite-amplitude effect. Longuet-Higgins 
concludes that the distribution of wave heights even ina 
Storm is well described by the Rayleigh distribution, pro- 
vided the rms amplitude, a, is estimated from the original 
record and not from the frequency spectrum as a = ( Zing * . 
The effect of finite bandwidth is estimated from a model 
assuming low background noise linearly superposed on a very 


narrow (delta function) spectrum. For narrow bandwidths, 
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he obtains the formula 


a*/2m, = 1 - 0.734 v* (2) 
where v is the rms spread of the noise about the mean fre- 
guency. Values of v2 corresponding to Pierson-Moskowitz 
(broad-band) spectra also give results in close agreement 
with observation. Therefore, the Rayleigh distributions 
calculated in this paper are parameterized using the rms 


wave height. 


Be LReNCATED PROBABILITY DISTRIBUTIONS 

In concept, waves are described by the joint distribu- 
tion of height, period (or equivalent wavelength) and direc- 
tion. To simplify the analysis, all authors assume a very 
narrow band frequency spectrum and a narrow directional spec- 
trum, so that all the wave heights of the distribution are 
associated with a single mean frequency and mean direction. 
Therefore, starting in deep water, the waves are described 
by the unaltered single-parameter Rayleigh distribution, 
with the implied assumptions. The deep water wave heights 
are transformed into shallow water waves using shoaling 
theory in which frictional dissipation is neglected. Even- 
tually the waves reach such shallow water that they start 
to break, with the largest waves breaking furthest offshore 
first. Wave breaking is simulated by truncating the tail 


of the Rayleigh distribution. 
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Peel linceoLStirloltion 
Collins (1970) was the first to apply the technique 
of a truncated distribution to describe the effects of wave 
breaking, using a Sharp cut-off with all broken waves equal 
to Hy which results in a delta function at Hh: Collins does 
not give an explicit formula for his distribution, but it 
would be the same as the described by Battjes (1974) (see 


below). He used linear shoaling theory and the breaking 


criterion after Le Mehauté and Koh (1967 


1/4 


= 1/7 2 
Hp/Hp = 0.76 tan’ 8(Hy/Lo) (3) 


where tan 8 is the bottom slope and Ho and Ly are the deep 


water wave height and length, respectively. For waves break- 


ing at an angle a, the bottom slope is actually tan 8 cos a 


ee 


b’ 


and H, should be replaced by H The various dis- 


0 0 0° 


tributions, and how they are truncated are showed schemati- 


COS a 


cally in Table I. 
Zeeebattjes Distribution 
Battjes (1974, 1978) again used a sharp cut-off of 
waves and applied the breaking criterion based on Miche's 
formula for the maximum height of periodic waves of constant 
f£OrmM , 


a) = 0.88/k tanh (y/0.88 kh) (4) 


where y is an adjustable coefficient. In shallow water (4) 


reduces to 
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Hy = yh. (5) 


He uses linear theory to shoal the waves. The probability 


density for breaking wave heights is given by: 


p (H) H/2H° exp[-1/2(H°/A")], for 0< H< H (6) 


p (H) exp[-1/2(H-/H*)16(H-H,), for H > 4H (7) 


b 


where H is the wave height parameterizing the truncated 
Rayleigh distribution by Battjes. All waves that have 

broken, or are breaking, assume the height prescribed by 
(7); this results in a delta function at the truncation 


height, H of the distribution (see Table I). 


ae 
Peewkue and Kuo Distribution 
Kuo and Kuo (1974) investigated the effect of break- 
ing on wave statistics using a conditional Rayleigh distri- 
bution sharply truncated, specified by the breaking wave 
height simply proportional to local water depth, equation 


(5). The conditional probability density function of wave 


heights, Py (A), is calculated using the following equation, 


2 pit) 
pee BH /0 <H <F ) = —te————— , for 0<H< H 
b b A 


J p(H)dH 
0 


b 


(8) 


= ag pee bomen > H 
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Describing the conditional probability in this manner re- 
Sees if the proportional redistribution of probability 
density associated with the broken or breaking waves over the 
range of H. This removes the delta function at the breaking 


wave height, H previously described by Collins and Battjes. 


b! 
Table I shows the original Rayleigh distribution with dotted 
lines, and the modified Rayleigh distributions in heavy lines 
after applying the cut-off to the tails using the breaking 
criterion. The distribution by Kuo and Kuo is more realis- 
tic but still results in a sharp cut-off of the distribution 
at the breaking heights. 
5. Goda Distribution 

Goda (1975) derived a more realistically truncated 
distribution, qualitatively anyway, by requiring a gradual 
cut-off of the distribution. He uses a shoaled Rayleigh 


distribution to describe the unbroken wave heights at shore- 


ward locations prior to applying his cut-off which is given 
by 

wm 
fe) 


_ 2 a 2 
po(H) = 4H/K¢H) exp (-2H*/K2H*) (9) 


where Ke is the shoaling coefficient. Goda (1975) calculated 
the wave shoaling using the nonlinear theory of Shuto (1974): 


which dictates the following: 


O < gHT*/h- < 30: Small Amplitude Theory (Op 


fl constant Ci 


30 < gHT“/h* < 50: Hh 
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50 < gHT*/h* < ©: Hh-/*[ VgHT“/h* - 2/3] = constant. (12) 


Goda assumes that wave breaking occurs in a range of wave 
heights between H. and Hos with varying probability. The 
probability density function of unbroken waves only is 


expressed as: 


p(X) = py fl) ; for X < X, (13) 
X-X, 

p(X) = P, (%) ~K) -X3 P(X) 7 fOr Xo ane = Ay (13) 

Pp, (%) = 0 stor Xy < X (15) 


where, normalized wave heights are defined by X = H/H 


7 


O 


X, = H,/H, and xX. = H/H The artifice of spreading breakers 


il 


over a range partly represents the inherent variability of 


‘oul 


breaker heights, and partly compensates for the simplifica- 
tion of using a single wave period in the estimation of 
breaker height. Broken waves generally have some height 
smaller than Xi - Since no theory is available for describing 
waves after they have broken, the heights of broken waves 

are assumed to be redistributed across the range and to be 
proportional to the unbroken waves as was done by Kuo and 

Kuo (1975). Therefore, the conditional probability density 


function for all heights is calculated by 


Pp, (X) 
I 


i) p, (X) dx 


zi 





where: 


xX eae 
uf -a X 


fp, (x dx 1- {1+a°x (x, -x,)}e +, (17) 
4 ; 


is a constant of proportionality applied to Pp, (%) [oanormnal— 
ize the pdf. In equation (17), the constant, a, is equal to 
V¥2/K_. 

The breaker height is estimated using the following 
formula, which is an approximate expression for Goda's 
breaker index (1970) based on laboratory data, 

L = H 
= A {1 -exp[-1.5 = -2(1+K tanP’g)]} , (18) 


H HL 
O O O 


es 
H 


where tan 8 denotes the bottom slope and the coefficients 
are assigned the following values for best-fitting to the 


index curves, 
A = 0.17, K = 15, and p = 4/3. 


The range of breaker height, X, -X is calculated 


i oy 


by assigning the following values for A: 


.- 2 = 
i et 2 cot XS 


The upper limit of A, was selected by considering the varia- 


bility of breaker heights, whereas the lower limit of A, was 


chosen simply as two-thirds of A The coefficients used 


1° 


ZZ 





are based on matching laboratory data taken on a 1/10 and 
1/50 beach slope and several field experiments. The Goda 
model is applied here to the experimental data described 
below. 
5. Summary 

The common idea of these studies is to cut-off the 
portion of wave height distribution beyond the breaker 
height, which is controlled by the water depth and other 
factors. The methods differ in the techniques of cut-off 


and the formulae used to define breaker heights. 
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Git. “EXPERIMENT 


Experiments were conducted at Torrey Pines Beach, San 
Diego, California, during November 1978, as part of the 
Nearshore Sediment Transport Study. At this site there is 
a gentle sloping, moderately sorted, fine-grained sandy beach. 
The beach profile shows no well-developed bar structure and 
is remarkably free from longshore topographic inhomogenei- 
ties. Winds during the experiments were light, and variable 
in direction. Shadowing by offshore islands and offshore 
refraction, limits the angles of wave incidence in 10m depth 
to less than 15°. During the experiments, Significant off- 
shore wave heights varied between 60 and 160 cm. The condi- 
tion of nearly normally incident, spilling (or mixed plung- 
ing~spilling) waves, breaking in a continuous way across the 


surf zone, prevailed during most of the experiments. 


A. INSTRUMENTS 

An extensive array of instruments was deployed to study 
nearshore wave dynamics. Measurements described here are 
from sensors located on an offshore transect from 10 m depth 
to across the surf zone (Fig. 1). The sensors were of three 
types: pressure (P), current (C) and surface~piercing-staff 
(W). 

The pressure sensors were Stathem temperature-compensated 
transducers with dynamic range of either 912-2316 g /om* or 


912~3720 g/em*. They were statically precalibrated and 
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postcalibrated by being lowered into a salt-water tank and 
were found guite linear; the gains differed by less than 2% 
between calibrations. 

Current meters were two-axis, Marsh-McBirney electro- 
magnetic, spherical (4 cm diameter) probes, with a three- 
pole output filter at 4 Hz. Precalibration and postcalibra- 
tion of current meters showed little change in replicate 
runs with steady or oscillating velocity fields. The uncer- 
tainty associated with using a single gain factor for all 
frequencies is roughly estimated at +5% in amplitudes (10% 
in variances). 

The wave staffs were dual resistance wires with low 
noise, high resolution, and good electronic stability. The 
accuracy of the wave staffs was about +3% based on repeata- 
bility of gain calibrations measured in the laboratory and 


in Siku. 


B. DATA ANALYSIS 

Sea surface elevation and wave velocity components were 
retrieved from sensors by telemetering to shore and there 
recorded on a special receiver/tape recorder, described in 
Mecail by Lowe et al., (1972). The sampling rate was 64 
samples/s which was reduced to 2 samples/s by digital low- 
pass filtering. Record lengths of approximately 68 minutes 
from each data set were analyzed. 

It was desired to examine only the sea-swell band of 


frequencies between 0.05 to 1.0 Hz (20 to 1 s periods). The 
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data were first linearly detrended to exclude effects of the 
rising and falling tides. The data were then high pass 
filtered with a cut-off frequency of 0.05 Hz (20 s period). 
The high pass filter used a Fast Fourier Transform algorithm 
to obtain the amplitude spectrum of the entire 68 minute 
record. The Fourier coefficients corresponding to 0 to 

0.05 Hz were used to synthesize a low frequency time series 
which was subtracted from the wave record. The limiting high 
frequency (Nyquist) was 1.0 Hz. 

The energy density spectra were calculated in a similar 
manner for wave and velocity measurements using the Fast 
Fourier Transform (FFT) algorithm. A cosine-squared taper 
data window was applied to the time series to minimize leakage. 

The highest maximum and lowest minimum of the surface 
elevation within a period interval defined, respectively, 
the crest and trough of a wave. A wave height H is defined 
as the total range of 7(t) in that interval, the time between 
two consecutive zero-up-crossings of c(t) (see Fig. 2). 

Since the average wave period was about 14 sec, the total 
number of waves in the 68 min. record was about 300, which 
gives reasonable wave statistics. The height statistics of 


mean wave height H, root mean square wave height Ht 


f 


mS 


Significant wave height Ay 3 (average of the heights of the 
1/3 highest waves), and Hy /10 (average of the heights of the 
1/10 highest waves), are calculated from the ordered set of 


wave heights. 
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The pdf and cumulative distributions of wave heights 
were calculated. The heights were normalized using the 
deep water root mean square value. Theoretical probability 
distributions were calculated using the Goda model and com- 
pared with measured distributions. 

A deep water reference wave height was calculated by 
measuring the energy using current meter C9 located at about 
4m depth and backing the energy out to deep water. Kinetic 
energy spectra were calculated from the measured horizontal 
velocity spectra, S,{f) and S(f); linear theory transfer 
functions were used to integrate the spectra over the water 


column, so that the average kinetic energy is 
_ 2 
<KE(£)> = fete ed) [S (£) +8 (£)] (ES?) 


where JH (£) |* is the transfer function that relates the 


velocity spectrum components to the kinetic energy, 


2 el: Sinh 2kh 


G2), | = fe) (20) 
aX cosh“k(h +z_) 


where Zo 1s the measurement elevation. Guza and Thornton 
(1980) showed, for these same experiments, that using linear 
theory transfer functions to calculate average kinetic energy 
gave reasonable results. To first order in energy, the 
average potential energy equals the average kinetic energy. 
The potential energy in deep water is obtained by applying 


linear shoaling transformation 


Fa 





7 2 
ge) ace, ee = [H. (£) | S10) 08) (21) 
where H_(f) = i Se is the linear shoaling coefficient. 


O 
The rms wave height is related to the PE and KE by 


<PE> = 


Co] F 
No 


5 ey ie ete =n | <KE(E)odr , (22) 
7 


from which the deep water rms wave height, denoted hence- 
becth by Hoe can be found. 

To take advantage of the large number of current meters, 
current data were used to infer wave heights. The velocity 
Signals were convolved using linear wave theory to obtain 
surface elevations. The complex Fourier spectra of the 
horizontal velocity components U(f), V(f) were first calcu- 
lated and vectorially added. The complex surface elevation 
spectrum, X(f), was calculated applying the linear wave 


theory transfer function, H(f) 


X(f) = H(£) -V(£) (23) 
where 
= sinh kh 
H(t) = a eosh kth Fz) Ai 


The complex surface elevation spectrum was then inverse 

transformed to obtain the surface elevation time series from 
which the wave height distribution is calculated. The entire 
68 min. record was convolved at one time in order to minimize 


the end effects which result in spectral leakage. 
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IV. RESULTS 


oee)6U LY PITCAL SPECTRA 

A broad range of wave and weather conditions were en- 
countered during the experiments. Typical velocity spectra 
for two days (Fig. 3, lower panel) include an example of 
very narrow band spectra calculated for November 20 for the 
current meters C22x and C23x which straddled the mean breaker 
line. In shallow water depths, the waves generally become 
more "peaky", resulting in increased energy at the harmonics. 
The presence of strong harmonics in the spectra indicates 
the importance of nonlinearities of the waves in shallow 
water. The spectral energy level decreases at all frequen- 
cies except at the very lowest from the deeper instrument 
C22x (heavy line) to the shallower instrument C23x due to 
breaking. The other typical spectra (Fig. 3, upper panel) 
are an example of combined sea and swell with a narrow band 
of energy at swell frequencies, but with broad band energy 


at higher frequencies. 


B. HEIGHT STATISTICS 

The wave height statistics for six days, inferred from 
wave staff and current meter data, are presented in Table II. 
The statistical parameters listed are: root-mean-square 


heights, H , Significant wave heights, Hy 3) average of 


rms 


the heights of the 1/10 highest waves, and the maximum 


Hy 107 


height, ax’ These parameters were obtained from 68 min. 
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records. The reference wave height He obtained by backing 
Piemeonernogy measured at current meter C9 out to deep water, 
the frequency peak of the spectrum and the depth for each 
instrument, are also presented. Most of the wave staff 
measurements were made inside the surf zone, after the waves 
have started breaking. 

Shoaling effects are observed in the data of Table II. 
The ee + Hy 7/3: and Hy io: increase as depth decreases, 
until they reach the breaking point after which the wave 
heights decrease as the depth decreases, due to breaking. 

The range of depths of the instruments is from about 
570 cm to 40 cm. The average peak frequency (peak frequency) 
varied little during the experiments and was about 0.07 sec + 
(Table II). The relative depth for the waves at the peak 
frequency is h/L < 1/25 in all cases so that they can be 


considered shallow water waves. 


C. COMPARISON OF EMPIRICAL WITH MODEL DISTRIBUTIONS 

The Goda model is compared with measured data. The model 
is first run using Goda's original coefficients. This model 
results in over-prediction of the Hms’ The coefficients 


are changed to optimize the model's description of the wave 


height distribution qualitatively, and H quantitatively. 


rms 


parameter was calculated from the second moment of 


The H 
rms 


the wave height distribution and is, therefore, a more sensi- 
tive parameter to describe the shape of the distribution 


(particularly the tail) than, say, first moments such as the 


mean H, So ey Ag 
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A comparison (Fig. 4) was made between the empirical 
wave height distribution normalized with the reference deep 
water wave height and Goda's model distribution for November 
20. Starting in deep water, it is observed that the model 
fits quite well the empirical distribution obtained from the 
pressure sensor P7. It is clearly seen that the wave heights 
in deep water are essentially Rayleigh distributed. [In 
shallow water the model over-predicts the tails and under- 
predicts the peaks. The smaller the depth the greater the 
errors (Figs. 4 and 5). Due to this over-prediction at the 
tails of the distribution the rms wave heights obtained by 
the model are larger than the measured ones. 

To obtain a better agreement between the measured and 
the calculated rms wave heights, the higher limit of breaking 
and the range of breaking of the model are changed. Goda 
(1975) calculated the shoaling using the nonlinear theory 
of Shuto (1974) (see equations 10, ll and 12). For all the 
ranges of wave heights, frequencies and depths used here, 
the values of ama /h- are calculated and in all cases analyzed, 
fall in the third category, equation(12). This law was used 
to calculate the nonlinear shoaling coefficient Ko: Goda 
assumed a range of breaking between Hl and H2 that takes 
into account the variability of breaker heights and the use 
of a Single frequency. The variable breaker height Hy isa 
function of the frequency, depth and bottom slope (see equation 
18). The values of the coefficients used by Goda in his 


breaking criterion are the following: Al = 0.18, A2 = 2/3 Al, 
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K = 15 and p = 4/3, empirically assigned to give the best- 
fit with observed wave heights. The coefficients used here 
to give a qualitatively better fit with the measured dis- 
tributions and quantitatively fit with eae are: Al = 0.136, 
A2 = 1/2 Al, K = 20; p was left equal to the previous value. 
The most sensitive coefficient is Al which was determined 
first. Values of Al were calculated for various distributions 
by first obtaining the higher breaker limit (Xl) from the 
measured empirical distributions as indicated by the maximum 
value of the distribution, and then calculating Al using equa- 
tion 18. The values of Al were then averaged to give a 
single representative value for all distributions. The 
coefficient A2 was chosen as 1/2Al, which results in a more 
symmetrical distribution as indicated by the results. The 
coefficient K, which weights the slope of the beach, was 
determined by trial-and-error testing of Goda's model for a 
variety of values looking for the best fit of all the distri- 
butions. The model is not very sensitive to changes in K. 
The comparison (Fig. 6) of the empirical distribution 
in deep water for November 20 with the distribution obtained 
applying the model used the new coefficients. In deep water, 
there is no notable difference from the original model's 
results (Fig. 4). In shallow water, the predicted values 
obtained with the modified model fit much better the 
empirical distributions than those obtained applying the 


original Gceda model. 
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As stated earlier, the coefficients Goda originally speci- 
fied were based on matching laboratory data for a two beach 
slopesof 1/10 and 1/50. The model was then applied to some 
field data, but unfortunately the beach geometry (slopes) 
were not given. The Torrey Pines Beach is approximately 
1/50 at the beach face and the wave climate is characterized 
by long period (~14 sec) swell. The reason for differences 
in the model coefficients needed to fit this data set com- 


pared with Goda's suggested coefficients is not known. 


D. COMPARISON OF RMS WAVE HEIGHTS 

The model calculated a values and measured values of 
a calculated directly from the wave heights are used to 
test how well the model works. Table III shows the measured, 
Rayleigh, and the calculated Goda and modified Goda root- 
mean-square wave heights corresponding to the wave staffs 
and current meters for six different days. 

Goda's model-predicted Hes values are, in general, 
larger than the measured ones (Fig. 9). In an attempt to 
explain the differences, _— values obtained from the 
model were plotted against depth, deep water wave height 
and Ursell number; but no meaningful correlation was obtained 
with any of these variables. 

Measured — and calculated values assuming the wave 
heights are Rayleigh distributed so that H,_. = 8m. were 


compared (Fig. 8a). Measured and assumed Rayleigh signifi- 


cant wave heights were also compared (Fig. 8b). The comparisons 
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show that the agreement between measured and Rayleigh sta- 
tistics are good for both small and large wave heights, 
inside and outSide the surf zone. The good comparisons with 
Rayleigh statistics suggests that the wave heights, although 
decreased by breaking, are still more nearly Rayleigh-dis- 
tributed than cut-off Rayleigh, as suggested by Goda. 

The change of the rms wave heights with depth, between 
the measured H ns and the calculated ones obtained by the 
model first applying nonlinear (heavy line) and then linear 
Shoaling (light line) were compared (Fig. 10). The model with 
nonlinear shoaling clearly Pe eaters the values, while 
the model with linear shoaling gives reaonable values, com- 
pared with measurements. The nonlinear shoaling "blows up" 
in very shallow water (<30 cm) and should be ignored. 

Measured and the calculated rms wave heights with depth 
obtained by applying linear and nonlinear shoaling using the 
modified coefficients were also compared (Fig. 10, lower 
panel). This figure illustrates that the modified Goda model 
fits better the data than the original does. The majority 
of the measured Bl fall between the two curves obtained with 
nonlinear and linear shoaling. Based on the choice of coeffi- 
cients, applying linear shoaling to Goda's model can give 
as good, or better, results as applying nonlinear shoaling. 

The rms wave height values, obtained using the modified 
model coefficients, were plotted against depth, deep water 


wave height and Ursell number; no obvious correlation could 
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be noted. It is found that the new predicted — values 
compared with measured ae have an error of less than +203 
at all depths (with the exception of one anomalous point) 


fic. 11). 


Peeeeevweve HEIGHT DISTRIBUTIONS USING CURRENT METERS 

As described earlier, the surface elevations were derived 
by linearly convolving the velocity records and wave height 
distributions calculated. The basis for applying this analy- 
Sis is the earlier work of Guza and Thornton (1980) where 
they showed, for this same data set, that linear theory spec- 
tral transformations could be used to calculate surface ele- 
vation standard deviations either from pressure meters or 
current meters with less than a 20% error, and typically less 
than 10%. Examples of the derived wave height distributions 
are considered for November 17 (Fig. 12). These measurements 
were made just outside the surf zone, at about the breaker 
point and inside the surf zone (current meters C23, C37 and 
C40 respectively). In general, the model overestimates the 
velocity derived wave height distributions more than the 
direct measurements. The reason for the discrepancy is that 
linear wave theory underestimates the surface elevations in 
convolving the velocities, particularly in the crest region 
of the waves. In other words, linear theory does not account 


for the finite amplitude of these highly nonlinear waves. 
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F. COMPARISON OF MODEL AND MEASURED CUMULATIVE DISTRIBUTIONS 

The cumulative exceedance of wave height distributions 
normalized with rms deep water wave height were calculated 
applying the modified Goda model (using nonlinear shoaling) 
and plotted with the measured cumulative distribution for 
comparison. The cumulative exceedance distribution empha- 
Sizes information in the tail of the distribution. For an 
example in shallow water, inside the surf zone (116 cm depth), 
there is a good agreement between the measured and predicted 
distributions with a slight underprediction by the model in 
the tail (Fig. 13). With sensors located just outside the 
surf zone, under-prediction of thetail is larger and over- 
prediction in the middle range occurs (Fig. 14). In general, 
there is a better agreement between the two distributions well 
inside the surf zone, e.g., wave staffs W38 and W4l, than 
those (W2l and W29) which were generally either at breaking 
or just outside the surf zone, the most nonlinear wave 


region. 
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V. CONCLUSIONS 


It is confirmed that the wave heights in deeper water 
(7 m) are Rayleigh distributed. 

Goda's model, using the empirical coefficients originally 
Suggested on the basis of laboratory and poorly specified 
field measurements, over-predicts the tail of the distribution 
and under-predicts the peaks. As a consequence, the predicted 
rms wave heights are larger than the measured. As the depth 
decreases, the errors of the predicted distributions increase. 
To get a better fit with the measured wave height distribu- 
tions, the coefficients in the breaking criterion used in 
Goda's model were modified. The values for Goda's breaking 
criterion giving the best fit to the measured data of this 


Study are: 
poe OeisG, eae = 72 Al and K = 20 


The percentage of error between the measured and the pre- 
dicted Has values from the model with these coefficients is 
les than +203. 

Linear shoaling was found to be as good as nonlinear 
Shoaling in applying Goda's model across the surf zone. 

Good comparisons were obtained between empirical and 
Rayleigh-derived statistics. This indicates that the wave 
heights, although decreased by breaking, are still more 
Rayleigh-distributed than the cut-off Rayleigh-distributed 


as suggested by Goda. 
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Goda's model was tested here with a large amount of 
field data for a variety of wave conditions on a 1/50 beach 
Slope. Further comparisons should be made for a variety of 


beach slopes and wave climates in order to test the general 


applicability of the model. 
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E(t) 


Figure 2. 





where x Ss a crest 


ms is a trough 


O is the zero-up crossing 


Definition sketch of zero-up-crossing 
wave heights. 
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Figure 3. Typical spectra measured during the experiments. 


44 





EMPIRICAL DISTRIGUTION CF WAVE MEIGHTS 
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Figure 4. Empirical distribution of wave heights compared 
with those predicted with Goda's model, starting 
in deep water (P7) and going into shallow water 
(W21, W38), 20 November 1978. 
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EMPIRICAL OISTAISUTION OF BiiVE HEIGHTS 
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Figure 5. Empirical distribution of wave heights compared 
with those predicted with Goda's model, starting 
in deep water (W21) and going into shallow water 
(W38, W411), 17 November 1978. 
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e 6. Empirical distribution of wave igh 
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model, 20 November 1978. 
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EMPIAICAL DISTHIGUTION CF WAVE HEIOHTS 
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Figure 7. Empirical distribution of wave heights compared 


with those predicted with the modified Goda's 
model, 17 November 1978. 
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MEASURED HAMS AT WAVE STAFFS VS RAYLEIGH HAMS 
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Figure 8a. Range of measured and Rayleigh 


Square wave heights. 


MEASURED At/3 AT WAVE STAFFS VS RAYLEIGH HI/3 
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Figure 8b. Range of meaSured and 
wave heights. 
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Figure 9. Correlation of measured rms wave heights 
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Figure 10. Comparison of the changes of Hymg with the Goda's 
model applying nonlinear (heavy line) and linear 
Sshoaling (light line). Upper figure illustrates 
Goda's original model; lower figure illustrates 
the use of the modified coefficients in applying 
Goda's model. 
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Figure ll. Percentage error of predicted (modified 
model) compared with measured rms wave 
heights. 
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Empirical distributions of wave heights 
obtained from current meters (C23, C37 and 
C40) compared with predicted wave heights 
calculated with the modified model, 

17 November 1978. 
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Figure 13. Comparison of measured cumulative exceedance 
distributions with predicted distributions 
(modified model), 20 November (W38) and 
17 November 1978 (W41). 
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Figure 14. Comparison of measured cumulative exceedance 
distributions with predicted distributions 
(modified model), 20 November (W21) and 
17 November 1978 (W38). 


oe 





BIBLIOGRAPHY 


Barber, N.F., "Ocean Waves and Swell," Lecture Published by 
Institution of Electrical Engineers, London, 22 pp., 
mo50. 


Battjes, J.A., “Computation of Set-Up, Longshore Currents, 
Run-UP and Overtoping Due to Wind-Generated Waves," 
1974. 


Battjes, J.A. and J.P.F.M. Janssen, "Energy Loss and Set- 
Up Due to Breaking of Random Waves," Proc. 16th Coastal 
enigameering Conf., Vol. I, pp. 569-587, 1978. 


Battjes, J.A., “Probabilistic Aspect of Ocean Waves," 


Report Now 77-2, Laboratory of Fluid Mechanics, Depart- 


ment of Civil Engineering, Delf University of Technology, 
1978. 


Cartwright, D.E. and M.S. Longuet~Higgins, "The Statistical 
Distributions of the Maxima of a Random Function," 


IgG. ROY. SOc., Ser A, 237, pp. 212-232, 1956. 


Collins, J.1I., “Long Shore Currents and Wave Statistics in 
the Surf Zone," Technical Report, Tetra Tech. No. 
meta 9—-2, February 1972. 


Collins, J.I., "Probabilities of Breaking Wave Characteris- 


mes, Proc. 12th Coastal Engineering Conf., pp. 399- 
Sine, 3970. 


G@@eeraoarty, S.K. and R.P. Cooley, "Statistical Distribu- 
tions of Periods and Heights of Ocean Waves," Jour. of 


Geophysical Research, Vol. 82, No. 9, pp. 1363-1368, 
Mamcho 20,1977. 


Forristal, G.Z., "On the Statistical Distribution of Wave 
Heights in a Storm," Jour. of Geophysical Research, 
eee NGew es, PD. 2353-c396, May’Z0, 19768. 


Goda, Y., “Numerical Experiments on Wave Statistics with 
Spectral Simulation," Report of the Port and Harbour 
Research Institute, Vol. 13, pp. 3-57, September 1970. 


Goda, Y., "A Synthesis of Breaker Indices," Jour. Civil 


Eng., Japan Soc. Civil Eng., No. 180, 1970 (in 


Japanese) 


56 





Goda, Y., "Irregular Wave Deformation in the Surf Zone," 


Seoccaleengineering jn Japan, Vol. 18, pp. 13-26, 1975. 


Guza, R.T. and E.B. Thornton, “Local and Shoaled Comparisons 
of Sea Surface Elevations, Pressures and Velocities," 


Moun. Of Geophysical Research, Vol. 85, No. C3, 
pp. 1524-1530, March 20, 1980. 


Kuo, C.T. and S.J. Kuo, “Effect of Wave Breaking on Statis- 
tical Distributions of Wave Heights," Proc. Conf. on 
Coastal Engineering, June 1974. 


LeMehauté, B. and R.C.Y. Koh, "On the Breaking of Waves 
Arriving at an Angle to the Shore," Journal of Hydraulic 
Research, Vol. 5, No. 1, 1967. 


Longuet-Higgins, M.S., "On the Distribution of the Heights 
of Sea Waves: Some Effects of Nonlinearity and Finite 


Bandwidth," Jour. of Geophysical Research, Vol. 85, 
feeeeeeo, pp. L5Sl9=-1523, March 20, 1980. 


Longuet-Higgins, M.S., "On the Statistical Distribution of 
the Heights of Sea Waves," Jour. of Mar. ReSearch, 
poe ll, No. 3, pp. 245-266, 1952. 


Longuet~Higgins, M.S., "On the Joint Distribution of the 
Periods and Amplitudes of Sea Waves," Jour. of Geo- 


physical Research, Vol. 80(18), pp. 2688-2694, 
20 June 1975. 


Longuet-Higgins, M.S., "The Effect of Nonlinearities on 
Statistic Distributions in the Theory of Sea Waves," 
Mouse Ot Fluid Mech., Vol. 17, pp. 459-480, 1963. 


Lowe, R.L., D.L. Inman, and B.M. Brush, "Simultaneous Data 
System for Instrumenting the Shelf," Proc. Conf. 


Pemdoeal Eng. 13th, pp. 95-ll2, 1972. 


Shuto, N., "Nonlinear Long Waves in a Channel of Variable 


Section," Coastal Engineering in Japan, Vol. 17, 
pee L=i2, 1974. 


Tayfun, M.A., "Linear Random Waves on Water of Nonuniform 


Pech, “CCean Bngineering Report No. 16, Part II, 


Department of Civil Engineering, University of Deleware, 
Newark, DE, August 1977. 


Tayfun, M.A., "Narrow-Band Nonlinear Sea Waves," Jour. of 
Geophysical Research, Vol. 85, No. C3, pp. 1548-1552, 
March 20, 1980. 


Dal 





magimecon, &.5., J.J. Galvin, F.L. Bub and P.D. Richardson, 
"Kinematics of Breaking Waves Within the Surf Zone," 


feoce oti Coastal Engineering Conf., pp. 461-476, 
woly 1976. 


Wood, W.L., "Wave Analysis System for the Breaker Zone," 


International Symposium on Ocean Wave Measurements and 
Analysis, Vol. I, pp. 774-789, 1974. 


58 





INITIAL DISTRIBUTION LIST 


Defense Technical Information Center 


Cameron Station 
Alexandria, VA 22314 


Library, Code 0142 
Naval Postgraduate School 
Monterey, CA 93940 


Chairman, Code 68 
Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93940 


Chairman, Code 63 
Department of Meteorology 
Naval Postgraduate School 
Monterey, CA 93940 


Professor E.B. Thornton, Code 68Tr 
Department of Oceanography 

Naval Postgraduate School 
Monterey, CA 93940 


Professor J. Wickham, Code 68Wk 
Department of Oceanography 
Naval Postgraduate School 
Monterey, CA 93940 


Galo Padilla Teran 
Teniente de Navio 


Instituto Oceanografico de la Armada 


Guayaquil, Ecuador 


Director 

Naval Oceanography Division 

Navy Observatory 

34th and Massachusetts Avenue NW 
Washington, D.C. 20390 


Commander 

Naval Oceanography Command 
NSTL Station 

payeeoe. bouls, MS 39522 


59 





10. 


i. 


2. 


ie 


14. 


>. 


IEG 


ne 


ie. 


19. 


Commanding Officer 

Naval Oceanographic Office 
NSTL Station 

Bay St. Louis, MS 39522 


Commanding Officer 
Fleet Numerical Oceanographic Center 
Monterey, CA 93940 


Commanding Officer 

Naval Ocean Research and Development Activity 
NSTL Station 

Bay St. Louis, MS 39522 


Office of Naval Research (Code 480) 

Naval Ocean Research and Development Activity 
NSTL Station 

meayeot. Louis, MS 39522 


Scientific Liaison Office 

Office of Naval Research 

Scripps Institution of Oceanography 
La Jolla, CA 92037 


Library 

Scripps Institution of Oceanography 
Eee, BOX 2367 

fa wolla, CA 92037 


Library 

Department of Oceanography 
University of Washington 
Seattle, WA 98105 


Library 

Crertsk 

P.O. Box 4803 

San YSidro, CA 92073 


Library 

School of Oceanography 
Oregon State University 
Corvallis, OR 97331 


Commander 

Oceanography Systems Pacific 
Box 1390 

Pearl Harbor, HI 96860 


60 





20 . 


ei. 


Ze. 


LS 


Chief, Ocean Services Division 

National Oceanic and Atmospheric Administration 
8060 Thirteenth Street 

eiiver Spring, MD 20910 


Director 
Instituto Oceanografico de la Armada 
Guayaquil, Ecuador 


Director de Educacion de la Armada 
Comandancia General de Marina 
Quito, Ecuador 


Director General del Personal 


Comandancia General de Marina 
MuLto, Ecuador 


61 




















193704 


Thesis 
me7os Teran 
e.. Weanstormation of 


waves across the surf 
zone. 


thesT2763 
ansto 


Ta cen ilSaiini ee | | 
IEA A 


DUDLEY KNOX LIBRARY 





